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ABSTRACT 

We  consider  the  problem  of  computing  the  shortest  path 
between  two  points  in  three-dimensional  space  bounded  by  a 
collection  of  convex  and  disjoint  polyhedral  obstacles.  Although 
this  problem  is  hard  to  solve  for  arbitrarily  many  obstacles,  it 
had  been  efficiently  solved  by  Mount  [Mo84]  (cf.  also  Sharir 
and  Schorr  [SS84])  for  a  single  convex  polyhedral  obstacle,  in 
time  C>(n^log  n),  where  n  is  the  number  of  faces  of  the  obstacle. 
•In  this  paper  we  generalize  this  technique  to  the  case  of  two 
convex  polyhedral  obstacles,  and  present  an  algorithm  which 
solves  this  problem  in  time  C>(n-'a(n)^("('')  ^log  n),  where  a(n) 
is  the  functional  inverse  of  Ackermann's  function,  and  is  thus 
extremely  slowly  growing.  This  result  is  achieved  by 
constructing  a  new  kind  of  Voronoi  diagram,  called  peeper's 
Voronoi  diagram,  which  we  introduce  and  analyze  in  this  paper. 


1.   Introduction 

Shortest  path  planning  in  Euclidean  space  bounded  by  obstacles  is  a 
problem  arising  in  robotics,  and  is  a  special  case  of  the  problem  of  planning 
optimal  collision  free  paths  of  a  given  robot  system  moving  amidst  a 
collection  of  known  obstacles,  between  an  initial  position  and  a  desired  final 
position.  In  general  this  problem  is  quite  hard,  and  the  notion  of  optimality  is 
not  easy  to  define,  but  if  the  robot  system  is  taken  to  be  a  single  moving 
point,  and  the  obstacles  are  assumed  to  be  a  collection  of  polyhedral  bodies, 
we  obtain  the  somewhat  easier  Euclidean  shortest  path  problem. 

In  the  2-dimensional  case,  the  shortest  path  between  two  points  X,  Y 
within  a  closed  plane  region  bounded  by  a  finite  collection  of  polygonal 
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obstacles,  is  a  polygonal  arc  whose  vertices  are  corners  of  the  given  polygonal 
obstacles,  and  it  can  be  easily  computed  by  discrete  graph  searching  in  time 
0(n^)  (cf.  [SS84],  [We84]). 

Suppose  that  V  is  a  closed  3-D  space  bounded  by  a  collection  of 
polyhedral  obstacles  {bd(y)  for  short),  and  let  X,Y  be  two  points  in  V.  Here 
the  shortest  path  from  X  to  Y  within  V  is  again  a  polygonal  path  whose 
corners  (except  for  X  and  Y)  lie  on  edges  of  bd{V)  (rather  than  at  corners  of 
bd{V)  as  in  the  2-D  case). 

To  reduce  the  problem  further  to  a  discrete  combinatorial  one,  Sharir 
and  Schorr  [SS84]  have  shown  that  shortest  paths  must  be  geodesic.  A  path  tt 
from  X  to  y  passing  through  the  lines  ^i.  .  .  .  ,?„  is  geodesic  if,  for  each 
1  =  1,  ...  ,n  the  path  tt  subtends  at  ^,  equal  incoming  and  outgoing  angles.  It 
is  shown  in  [SS84]  that  there  exists  a  unique  geodesic  path,  hence  a  unique 
shortest  path  -rr  from  X  to  7  which  is  constrained  to  pass  through  a  given 
sequence  of  lines  ^i,  .  .  .  ,^„. 

This  property  enables  us  to  consider  the  3-D  path  finding  problem  as 
essentially  solved,  if  we  can  specify  the  sequence  of  obstacle  edges  and 
corners  through  which  the  desiieid  shortest  path  must  pass.  In  this  setting  the 
problem  is  reduced  to  a  purely  combinatorial  one,  and  can  be  solved  in  0{n") 
steps,  each  step  consisting  of  finding  a  shortest  path  constrained  to  pass 
through  some  specific  sequence  of  ol'stacie  edges  and  comers. 

Computing  the  points  of  contacts  (and  thus  also  the  length)  of  a  geodesic 
path  passing  through  a  given  sequence  of  edges  also  seems  to  be  an 
intractable  problem.  Sharir  and  Schorr  show  that  the  exact  solution  of  this 
problem  involves  solution  of  a  system  of  high-degree  algebraic  equations, 
mainly  because  any  two  adjacent  edges  ^/_i  and  ^i  in  the  sequence 
^=(^1,  .  .  .  ,^„)  of  wall  edges  in  bd(y)  through  which  tt  passes,  can  be  skew 
to  one  another.  Sharir  and  Schorr  show  that  this  system  of  equations  is 
solvable  in  time  doubly  exponential  in  the  number  of  edges,  a  result  which 
has  recently  been  improved  to  a  single  exponential  method  by  Reif  and  Storer 
[RS85]. 

1.1.   Shortest  Paths  Along  a  Convex  Polyhedron 

In  contrast  to  the  difficulties  in  solving  the  general  3-D  shortest  path 
problem,  the  special  case  involving  a  single  obstacle  which  is  a  convex 
polyhedron  is  shown  in  [SS84]  and  [Mo84]  to  be  solvable  in  polynomial  time. 
We  briefly  review  the  basic  properties  of  shortest  paths  in  this  case,  and 
outline  Mount's  procedure  for  their  calculation.  First  note  that  the  shortest 
path  between  two  points  X,y  on  the  surface  5  of  a  convex  polyhedron  K, 
which  avoids  the  interior  of  K,  Hes  entirely  along  S.  Thus,  in  this  case,  any 
two  adjacent  edges  ^/_i  and  ^i  which  are  crossed  by  a  shortest  path  it  along 
5,  cannot  be  skew  to  one  another,  thereby  making  the  problem  of  computing 
points  of  contact  of  tt  with  the  edges  of  K  immediately  solvable  as  follows. 

Definition: 

(a)    A  point  Z  €  5  is  called  a  "ridge  point"  (relative  to  the  fixed  source  point 
X)  if  there  exist  at  least  two  shortest  paths  from  X  to  Z  along  5.  We 
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denote  by  R  the  set  of  all  ridge  points  in  S. 

(b)  For  each  point  Z  €  S-R,  let  -niZ)  denote  the  unique  shortest  path  from 
XtoZ  along  5. 

Let  Z  ^  S  and  suppose  that  the  sequence  C=(Ci.  .  .  .  ,C„)  of  edges  of  K 
through  which  'n{Z)  passes  has  already  been  found.  Then  we  also  know  the 
corresponding  sequence /o,  •••./„  of  faces  of  K  through  which  •7t(Z)  passes, 
where//  contains  the  two  edges  \t  and  ^/+i,  for  /=1,  .  .  .  ,n— 1  and  where /q 
contains  X  and  ^i  and  /„  contains  ^„  and  Z. 

We  then  unfold  the  collection  of  faces  /q,  .  .  .  ,/„  so  as  to  make  them  all 
lie  in  the  same  plane  L  (without  overlapping  one  another).  The  required  path 
ir(Z)  unfolded  to  L,  is  then  simply  the  straight  segment  XZ.  The  points  of 
intersection  of  this  segment  with  the  unfolded  edges  4i,  .  .  .  ,C„  are  then 
readily  determined,  and  can  be  easily  transformed  back  to  the  original 
polyhedron. 

The  main  issue  still  left  is  the  combinatorial  problem  of  finding  the 
sequence  of  edges  through  which  the  desired  shortest  path  must  pass.  Mount 
[Mo84]  tackles  the  problem  in  which  the  source  point  is  fixed  and  the 
destination  point  is  arbitrary,  by  defining  an  equivalence  relation  on  the 
(non-ridge)  points  of  the  polyhedron,  making  two  such  points  equivalent  if 
the  shortest  paths  to  these  points  traverse  the  same  sequence  of  faces.  An 
equivalence  class  of  this  relation  forms  a  region  on  a  face  of  the  polyhedron 
called  a  slice.  The  slice  boundaries  consist  of  edges  of  the  polyhedron  and  of 
ridges.  .".:.„:>  • 

Mount  presents  an  algorithm  which  partitions  the  faces  of  the 
polyhedron  into  slices,  in  time  0{n^\og  n)  .  The  algorithm  is  based  on  the 
observation  that  all  slices  arise  as  Voronoi  cells  of  certain  sets  of  points,  and 
will  be  briefly  reviewed  later  on. 

In  this  paper  we  study  a  generalized  version  of  the  3-D  shortest  path 
problem  in  the  presence  of  just  two  disjoint  convex  obstacles. 

The  main  part  of  this  paper  treats  the  problem  of  computing  the  shortest 
path  between  two  points  lying  respectively  on  two  convex  polyhedra  -  a 
source  polyhedron  and  a  destination  polyhedron.  We  present  a  shortest  path 
algorithm  that  runs  in  0(/z^a(n)'^^°('')^Mogn)  time,  where  a{n)  is  the 
functional  inverse  of  Ackermann's  function,  and  is  thus  extremely  slowly 
growing. 

Like  Mount's  algorithm  our  algorithm  starts  by  partitioning  the  faces  of 
each  polyhedron  into  slices.  We  then  make  use  of  the  fact  that  the  required 
shortest  path  must  be  the  concatenation  of  subpaths  having  the  following 
form:  - 

(a)  a  shortest  subpath  lying  entirely  on  the  source  polyhedron's  surface. 

(b)  a  straight  line  segment  connecting  a  take-off  point  on  the  source 
polyhedron's  surface  with  a  landing  point  on  the  target  polyhedron's 
surface. 

(c)  a  shortest  subpath  lying  entirely  on  the  target  polyhedron's  surface. 
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The  basic  idea  behind  our  algorithm  is  that  the  take-off  point  and  the 
landing  point  fully  determine  the  required  shortest  path,  because  the  problem 
can  then  be  split  into  two  already  solved  subproblems,  each  of  which  deals 
with  calculation  of  a  shortest  path  along  a  single  convex  polyhedron.  Thus 
the  problem  is  reduced  to  finding  those  two  crucial  points. 

The  paper  is  organized  as  follows.  In  Section  2  we  review  the  basic 
properties  of  shortest  paths  along  a  single  convex  polyhedron,  and  briefly 
describe  Mount's  technique  for  the  calculation  of  such  paths.  In  section  3  we 
analyze  the  shortest-path  problem  involving  two  convex  obstacles,  and 
provide  an  initial  crude  solution  to  this  problem.  This  solution  is  refined  and 
improved  in  Section  4,  by  the  use  of  a  certain  modified  3-D  Voronoi  diagram, 
called  peeper's  Voronoi  diagram,  which  is  introduced  and  analyzed  in  this 
section.  The  analysis  of  this  diagram  involves  Davenport-Schinzel  sequences, 
and  uses  recent  results  on  the  maximal  length  of  such  sequences.  We  then  use 
peeper's  Voronoi  diagrams  for  solving  our  shortest-path  problem  in  time 
C>(n^a(n)^(°'(")  ^log  n),  which  is  almost  proportional  to  0{nHog  n).  In  section 
5  we  generalize  the  problem  to  cases  where  the  source  and  destination  points 
of  the  shortest  path  lie  anywhere  in  the  closed  three  dimensional  region 
bounded  by  the  obstacles. 

2.   Mount's  Algorithm  for  Shortest  Paths  Along  a  3-D  Convex  Polyhedron 

2.1.   The  Structure  of  Shortest  Paths 

Consider  a  convex  polyhedron  K  in  3-space  with  n  vertices.  By  Euler's 
formula,  the  number  of  edges  and  faces  of  the  polyhedron  is  0{n).  Let  sq 
denote  a  fixed  source  point  from  which  we  wish  to  calculate  shortest  paths  to 
any  other  point  on  the  surface  S  of  K,  and  let/o,  •  .  •  ,fm  denote  the  faces  of 
S  where /o  is  the  face  containing  sq. 

Definition:  hctfijj  be  adjacent  faces  on  S.  The  planar  unfolding  offj  relative 
to  fi  is  an  orthogonal  transformation  which  maps  the  points  of  fj  into  the 
plane  of /y  by  unfolding  fj  about  the  edge  between  these  faces,  into  the  plane 
containing//  (so  that  it  does  not  overlap  with/).  Similarly,  if /i,  .  .  .  ,/,  is  a 
sequence  of  faces  such  that  /  and  fi+i  are  adjacent  in  S  for  each 
i  =  l,  .  .  .  ,r— 1,  then  the  planar  unfolding  of  this  sequence  relative  to  fi  is 
obtained  by  unfolding  these  faces  one  at  a  time,  about  the  edges  which 
separate  them,  until  they  all  lie  in  the  plane  containing  /^  (with  no  'two 
adjacent  faces  overlapping  one  another). 

We  recall  the  following  properties  (for  which  see  [SS84],[Mo84]): 

(1)  Suppose  we  know  the  sequence  of  edges  crossed  by  the  shortest  path 
from  Sq  to  some  point  Y,  then  that  shortest  path  can  be  determined  by 
taking  the  planar  unfolding  of  the  corresponding  sequence  of  faces 
incident  to  these  edges,  relative  to  /o,  then  by  drawing  the  straight 
segment  connecting  sq  to  the  (unfolded)  point  Y,  and  finally  by  folding 
this  straight  segment  back  into  5. 

(2)  Two  shortest  paths,  emanating  from  the  common  source  ^o  do  not  cross 
each  other. 
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(3)   A  shortest  path  cannot  pass  through  a  vertex  of  the  polyhedron. 

Definition:  A  point  p  is  a  source  image  for  a  point  x  on  some  face  /  if  p  lies 
in  the  plane  containing  /  and  if  the  segment  px  is  the  planar  unfolding 
(relative  to  /)  of  a  shortest  path  from  the  source  point  sq  to  j:  . 

Planar  unfolding  of  the  shortest  path  from  ^o  to  a  point  x,  can  also  be 
done  symmetrically  in  the  plane  containing  /q.  Doing  so  for  every  point  on 
the  polyhedron's  surface,  yields  an  unfolding  of  S  into  the  plane  containing  /o 
so  that  shortest  paths  become  straight  rays  emanating  from  the  source  sq  (see 
Fig.  2.1).   Note  that  in  this  unfolding  faces  will  be  split  up  along  the  ridges. 
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Figure  2.1  Planar  Layout  of  a  Polyhedron 

Definition:  A  slice  is  a  subset  of  S  —  R  consisting  of  all  points  Z  such  that 
the  shortest  paths  'n{Z)  all  traverse  the  same  sequence  of  edges  of  K. 

Note  that  the  boundaries  between  slices  are  formed  by  the  edges  of  the 
polyhedron  and  by  the  ridge  points. 

Lemma  2.1  ([SS84]) 

(1)  The  set  of  ridge  points  is  the  union  of  finitely  many  line  segments  called 
ridges. 

(2)  The  ridges  form  a  continuous  tree  on  the  polyhedron  whose  leaves  are 
vertices  of  the  polyhedron  and  whose  vertices  of  degree  2  are 
intersections  of  ridges  and  polyhedron  edges. 

(3)  Each  slice  is  a  convex  polygon. 

(4)  Shortest  paths  do  not  pass  through  ridges. 
Lemma  2.2  ([Mo84]) 

(1)  There  are  at  most  n  slices  on  a  given  face. 

(2)  There  are  0{n^)  slices  over  the  entire  polyhedron. 

(3)  There  are  0{n^)  ridges  over  the  entire  polyhedron. 


-6- 

(4)  For  all  points  in  the  same  slice,  the  planar  unfolding  of  the  shortest  paths 
to  them,  relative  to  the  face  containing  the  slice,  are  straight  lines 
emanating  from  a  single  point,  which  is  an  unfolded  image  of  the  source. 

Having  calculated  these  slices.  Mount  transforms  the  shortest  path 
problem  into  a  standard  point  location  problem.  Once  we  know  which  slice 
contains  a  given  target  point  x,  the  shortest  path  to  x  can  be  easily 
determined  by  drawing  a  straight  segment  from  x  to  the  corresponding 
unfolded  source  image,  and  by  folding  this  segment  back  over  the 
polyhedron. 

2.2.   Computation  of  the  Slices 

Let  /  be  a  face  of  K.  The  preceding  results  imply  (cf.  [SS84],  [Mo84]) 
that  there  exist  at  most  n  source  images  pi,  .  .  .  ,/?„  in  the  plane  containing/ 
such  that  the  slices  within  /  are  precisely  the  intersections  of  /  with  the  cells 
of  the  standard  (EucUdean)  Voronoi  diagram  of  the  images  pi,  .  .  .  ,p„. 

This  implies  that  if  e  is  an  edge  of  K,  then  the  ridges  split  e  into  at  most 
n  subintervals,  at  points  which  are  the  intersections  of  e  with  the  Voronoi 
diagram  induced  by  all  source  images  that  correspond  to  shortest  paths  which 
reach  points  on  e.  We  therefore  call  each  such  subinterval  /  a  Voronoi 
subinterval,  and  denote  it  by  Vor{p  ,e),  where  p  is  the  source  image  such  that 
all  the  (vmf olded)  shortest  paths  to  points  on  /  reach  /  from  p . 

Mount's  algorithm  simply  constructs  all  these  restricted  one-dimensional 
Voronoi  diagrams  along  edges  of  K  incrementally,  by  adding  source  images 
to  these  diagrams  one  at  a  time.  Whenever  the  algorithm  adds  a  new  source 
image  to  the  Voronoi  diagram  on  some  edge  e,  it  unfolds  and  propagates  this 
image  to  the  other  edges  of  K  that  are  co-planar  with  e,  as  a  potential  new 
source  image  for  these  edges.  The  algorithm  also  maintains,  for  each 
candidate  source  image  p'  associated  with  an  edge  e' ,  the  length  of  the 
smallest  shortest  path  that  can  reach  e'  from  p'  (and  not  be  beaten  by 
shortest  paths  from  other  source  images  already  found  for  e').  These  lengths 
are  maintained  in  a  priority  queue,  so  that  the  next  incremental  step  of  the 
algorithm  always  chooses  the  candidate  source  image  with  the  smallest 
associated  path  length. 

Mount  shows  that  each  step  in  this  algorithm  can  be  accomplished  in 
essentially  (9(log  n)  time,  and  since  the  total  number  of  slices  is  O(n^),  the 
algorithm  calculates  all  Voronoi  subintervals  on  edges  of  K  in  total  time 
0{n^\og  n).  Once  these  intervals,  and  their  corresponding  source  images, 
have  been  determined,  the  slices  within  each  face  f  of  K  can  then  be 
calculated  from  the  Voronoi  diagram  of  the  source  images  "reaching"  /.  Each 
such  diagram  can  be  calculated  in  time  0(n  log  n),  so  that  the  entire 
collection  of  slices  can  also  be  calculated  in  time  0{n^\og  n).  For  more  details 
concerning  this  procedure,  see  [Mo84]. 

Summing  up,  we  have 
Theorem  2.3  ([Mo84]) 
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Given  a  convex  polyhedron  with  n  vertices  and  a  source  point  sq  on  it,  the 
shortest  path  from  sq  to  any  query  point  on  the  surface  of  the  polyhedron  can 
be  determined  in  0{k  +  \ogn)  time  after  C>(n^log  n)  preprocessing  from  a 
structure  requiring  0(n^)  storage.  Here  k  is  the  number  of  faces  traversed  by 
the  required  shortest  path. 

Remark:  By  using  more  sophisticated  data-structures,  Mount  was  able  to 
further  reduce  the  space  requirement  of  the  algorithm  to  0{n  log  n). 
However,  in  what  follows  we  will  use  the  simpler  algorithm  as  outlined 
above. 

3.  Shortest  Paths  Between  Two  Convex  Polyhedra 

Let  A,  5  be  two  disjoint  compact  convex  polyhedra  in  3-space.  We  will 
denote  by  m  and  n  the  number  of  edges  in  A  and  B  respectively,  and  assume 
without  loss  of  generality  that  n>m. 

Let  Sq  denote  a  source  point  lying  on  the  surface  of  polyhedron  A,  and 
let  Z  be  a  destination  point  lying  on  the  surface  of  polyhedron  B.  The 
problem  studied  in  this  section  is  that  of  finding  the  shortest  path  from  sq  to 
Z  which  avoids  the  interiors  of  both  A  and  B.  We  call  A  the  source 
polyhedron  and  B  -  the  destination  polyhedron. 

For  simplicity  of  presentation,  we  assume  that  both  points  sq  and  Z  lie  in 
the  interior  of  a  face.  This  assumption  can  always  be  enforced  by  making  an 
infinitesimal  perturbation  in  the  location  of  these  two  points  if  they  happen  to 
lie  on  an  edge  or  at  a  vertex. 

In  contrast  to  the  case  of  a  single  polyhedron,  we  will  only  consider  here 
the  case  in  which  the  source  and  the  destination  points  are  both  fixed.  This  is 
because  the  more  complex  geometry  of  shortest  paths  in  the  case  of  two 
polyhedral  obstacles  makes  it  much  more  expensive  to  pre-process  the  data 
before  knowing  the  destination  point  (as  is  done  e.g.  in  [Mo84]). 

3.1.  The  Structure  of  Shortest  Paths 

Lemma  3.1:  The  shortest  path  tt  from  sq  to  Z  is  the  concatenation  of  the 
following  three  subpaths: 

(a)  A  shortest  subpath  along  the  surface  of  A  from  sq  to  a  "take-off"  point  t. 

(b)  A  straight  line  segment  in  empty  space,  connecting  the  take-off  point  t  to 
a  "landing"  point  /  on  the  surface  of  B . 

(c)  a  shortest  subpath  along  the  surface  of  B  from  the  landing  point  /  to  the 
destination  Z. 

Proof:  By  standard  arguments  (cf.  fSS84])  it  is  a  polygonal  path  each  of 
whose  comers  (except  sq  and  Z)  lies  on  an  edge  of  A  or  on  an  edge  of  B . 
Furthermore,  since  A  and  B  are  convex,  any  subpath  of  tt  between  two  points 
on  A  must  lie  entirely  along  the  surface  of  A,  and  similarly  for  B.  These 
observations  plainly  imply  the  lemma.  □ 

Remark:  In  cases  where  the  take-off  point  coincides  with  ^o,  the  first  subpath 
degenerates  to  a  single  point;  similarly  if  the  landing  point  coincides  with  Z, 
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the 


third  subpath  degenerates  to  a  single  point.  Since  A  and  B  are  disjoint 
polyhedra,  the  second  segment  can  never  degenerate  to  a  point. 

Lemma  3.1  implies  that  once  we  have  found  the  take-off  and  the  landing 
points,  we  can  calculate  the  required  shortest  path  by  separately  finding  the 
shortest  path  along  the  surface  of  A  from  the  source  point  j'q  to  r,  and  the 
shortest  path  along  the  surface  of  B  from  the  "source"  point  Z  to  /.  Since  A 
and  B  are  convex  polyhedra,  each  of  these  two  subpaths  can  be  calculated 
using  Mount's  algorithm  for  a  single  convex  polyhedron. 

3.2.   Finding  the  Take-off  and  Landing  Points 

Lemma  3.1  restricts  the  search  for  the  take-off  and  landing  points  to  the 
edges  of  A  and  B  plus  the  source  point  sq  and  destination  point  Z.  To  reduce 
the  problem  further  to  a  purely  discrete  combinatorial  one,  we  proceed  as 
follows. 

Let  ^1,  ...  ,^/  be  a  sequence  of  edges  of  A,  traversed  by  some  shortest 
path  along  the  surface  of  A  from  sq  to  some  point  in  the  interior  of  ^,. 
Similarly,  let  t^i,  .  .  .  ,T|y  be  a  sequence  of  edges  of  B,  traversed  by  some 
shortest  path  along  the  surface  of  B  from  Z  to  some  point  in  the  interior  of 
Ti^.  We  denote  (^i,  .  .  .  ,^,)  by  ^  and  (tii,  .  .  .  ,-x\j)  by  t]. 

Let  I,  denote  the  sequence  (^i,  .  .  .  ,^i,'r\j,'t\j-i,  .  .  .  ,r[i).  Our  current 
goal  is  to  find  out  whether  I,  can  be  a  candidate  for  the  sequence  of  edges 
traversed  by  the  desired  shortest  path.  If  so  then  we  would  also  like  to 
compute  the  corresponding  take-off  and  landing  points  t  and  /,  in  order  to 
determine  the  shortest  path  to  Z  constrained  to  traverse  4-  We  denote  this 
path  by  t:^. 

Note  that  on  every  path  traversing  ^,  the  take-off  point  lies  on  ^,  and  the 
landing  point  lies  on  T|y.  Also  by  the  results  reviewed  above,  the 
corresponding  take-off  and  landing  points  are  uniquely  determined  by  I,. 

Lemma  3.2:  The  shortest  path  ir  cannot  pass  through  any  vertex  of  the 
polyhedra  A,  B. 

Proof:  Qearly,  the  first  and  third  subpaths  of  tt,  as  in  Lemma  3.1,  cannot 
pass  through  a  vertex  of  A  and  B  respectively.  We  also  claim  that  the  take- 
off and  landing  points  themselves  cannot  be  vertices.  Suppose  to  the  contrary 
that  the  take-off  point  t  does  coincide  with  some  vertex  of  A.  Let  Hull  (A, I) 
denote  the  convex  hull  of  the  union  of  polyhedron  A  with  the  landing  point  /. 
Assuming  that  the  points  sq,1  are  not  visible  from  each  other  (or  else  t  would 
have  coincided  with  ^o)»  the  shortest  path  from  .yo  to  /  along  the  surface  of 
Hull(A,l)  must  be  the  concatenation  of  the  first  two  subpaths  of  tt, 
contradicting  the  fact  that  a  shortest  path  along  a  single  convex  polyhedron 
cannot  pass  through  a  vertex.   The  proof  for  /  is  completely  symmetric,  o 

We  next  transform  the  subproblem  of  finding  a  geodesic  path  tt^  passing 
through  4  into  the  context  of  Mount's  algorithm.  Suppose  the  edges  it-\,^t 
lie  on  a  face  /„  of  A,  and  the  edges  'r]j-i,T[j  lie  on  a  face  f^  of  B.  The 
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sequence  ^  induces  a  Voronoi  interval  Vor(p  ,^t)  or*  ?/.  whose  associated 
source  image  p  is  the  image  of  sq  under  the  planar  unfolding,  relative  to  /„, 
of  all  faces  of  A  sharing  the  edges  in  ^.)  Similarly,  the  sequence  r\  induces  an 
interval  Vor(q,r\j)  on  T|y,  where  q  denotes  the  image  of  the  destination  Z 
resulting  from  the  planar  unfolding,  relative  to/^,  of  the  faces  of  B  sharing 
the  edges  in  t|  . 

Assume  for  the  moment  that  each  pair  of  points  in  Vor(p,^i)  and  in 
Vor(q,r\j)  are  visible  from  each  other.  We  then  seek  two  points  t,l  lying 
respectively  in  these  two  intervals,  which  minimize  the  total  length  of  the 
path  pt  W  tl  W  Iq.  However,  since  this  path  consists  of  three  segments,  its 
corners  can  be  calculated  in  constant  time  using  the  technique  described  in 
[SS84].  Note  that  this  technique  yields  points  r,  /  which  can  he  anywhere  on 
the  lines  containing  ^/  and  iiy  respectively.  Nevertheless  we  must  have  the 
following  obvious 

Lemma  3.3:  Let  t,l  denote  the  take-off  and  landing  points  in  tt^.  Then  a 
necessary  condition  for  ir^  to  be  the  required  shortest  path  tt  is  that 
t  €  Vor(p,^/)  and  /  C  Vor{q,'^j)  (and  neither  t  nor  /  are  vertices). 

Proof:  If  r  £  Vorip  Xi)  then  t  is  not  contained  in  the  slice  induced  by  ^  on  the 
face  /„  of  A.  Hence  the  shortest  path  from  sq  to  t  does  not  pass  through  ^. 
Similarly,  if  /  £  'Vor{q,'r\j)  then  the  shortest  path  from  Z  to  I  cannot  pass 
through  T| .  □ 

Hence,  after  calculating  the  points  r,/  we  check  that  the  conditions  in 
Lemma  3.3  are  satisfied.  Only  then  ttj  is  considered  as  a  feasible  candidate 

for  TT. 

Ignoring  possible  non-visibihty  of  edges  of  A  from  edges  of  B,  the 
preceding  analysis  thus  suggests  the  following  crude  and  straightforward 
technique  for  finding  the  shortest  path  tt.  Since  each  edge  of  A  or  of  5 
contains  at  most  0{n)  Voronoi  subintervals,  there  are  at  most  0(/i'*)  pairs  of 
such  Voronoi  intervals,  one  on  an  edge  of  A  and  the  other  on  an  edge  of  B. 
For  each  of  these  pairs  we  check  whether  the  shortest  path  traversing  the 
corresponding  concatenated  sequence  of  edges  of  A  and  of  B  is  a  candidate 
for  IT  (in  the  sense  of  the  necessary  condition  in  Lemma  3.3).  If  so,  then  we 
compute  its  length  in  0(1)  time,  tt  is  then  the  shortest  path  among  all  those 
0{n^)  candidate  paths. 

This  crude  algorithm,  which  runs  in  time  0{n'^),  will  be  improved  later 
on  in  Section  4.  But  before  doing  this,  we  first  have  to  analyze  the  visibility 
between  edges  of  A  and  edges  of  5,  which  we  have  so  far  been  ignoring. 

3.3.  Visibility  Between  Two  Edges  on  Two  Convex  Polyhedra 
Definitions: 

(a)  Two  points  x^A,y^B  are  visible  from  one  another  if  the  segment  xy  does 
not  intersect  the  interiors  of  A  and  of  B . 

(b)  A  line  segment  /  is  obscure  to  a  point  x ,  if  there  is  no  interior  point  of  / 
visible  from  j:  . 
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(c)  A  line  segment  /  is  visible  from  point  x,  if  every  point  along  /  is  visible 
fromx. 

(d)  Two  line  segments  /i,/2  are  obscure  to  each  other,  if  there  are  no  two 
interior  points  of  l^  and  I2  visible  from  each  other. 

(e)  Two  line  segments  1^,12  are  visible  from  each  other,  if  every  point  along 
/i  is  visible  from  every  point  along  I2. 

Lemma  3.4:  Let  a;  be  a  point  on  A  lying  in  the  interior  of  a  face  f^  and  iet  e 
be  some  edge  of  B.  Then  e  consists  of  at  most  two  subintervals,  one  visible 
from  X  and  the  other  obscure  to  x. 

Proof:  LetfiJj  denote  the  faces  of  B  incident  to  e,  and  let  UiJlj  be  the  two 
planes  containing  fiJj  respectively.  By  convexity  of  fi,  11/  (resp.  IL)  is  a 
supporting  plane  of  B.  Let  Hf  (resp.  Hj)  denote  the  closed  half  space 
determined  by  11,  (resp.  IL)  whose  interior  is  disjoint  from  B,  and  let  Hf 
(resp.  Hj)  denote  the  complementary  open  half  space. 

Qearly  e  is  obscure  to  every  point  in  HfHHj,  and  is  visible  from  every 
point  in  HfUHJ,  disregarding  obstacles  other  than  B  (see  Fig.  3.1).  Theiefore 
if  X  lies  in  HfOHJ,  e  is  obscure  to  x  and  we  are  done. 

Otherwise  let  H^  denote  the  plane  containing  /t .  By  convexity  of  A ,  11^  is 
a  supporting  plane  of  A .  K  IIj^  does  not  intersect  e ,  then  either  e  lies  wholly 
in  the  half-space  bounded  by  Iljfc  and  containing  A,  in  which  case  e  is  again 
obscure  to  x,  or  else  e  lies  wholly  in  the  other  side  of  ITjt,  in  which  case  it  is 
clearly  visible  from  x.  Finally  if  fl^  intersects  e,  it  splits  e  into  two  portions, 
one  obscure  to  x  and  the  other  visible  from  x .  □ 


Figure  3 . 1 
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Lemma  3.5:  Let  g,e  be  two  edges  of  A,B  respectively,  and  let  x  be  a  point 
lying  on  g.  Then  either  e  is  obscure  to  a:,  or  else  there  exists  a  (possibly 
empty)  subsegment  e»  of  e  obscure  to  x,  whose  end  points  depend  on  e  and 
g,  but  not  on  the  choice  of  x^g.  Thus  Cg  splits  e  into  at  most  two  remaining 
portions,  each  of  which  is  visible  from  every  point  on  g  to  which  e  is  not 
(totally)  obscure. 

Proof:  Let  /t/,  (resp,  fiJj)  denote  the  faces  of  A  (resp.  B)  incident  to  g 
(resp.  e),  and  let  11^,11/  (resp.  TliJlj)  be  the  two  planes  containing /t/;  (resp. 
fiJj)  .  By  convexity  of  A,  rijt  (resp.  FI/)  is  a  supporting  plane  of  A.  Let  Hi- 
(resp.  Hf)  denote  the  closed  half  space  bounded  by  Fljt  (resp.  IT/)  whose 
interior  is  disjoint  from  A,  and  let  H^  (resp,  Hf)  denote  the  complementary 
open  space. 

Suppose  for  the  moment  that  the  faces  of  B  were  transparent,  so  that  e 
could  be  seen  behind  them.  Three  cases  can  arise: 

(a)  e  lies  entirely  in  H^dHf.  In  this  case  e  is  obscure  to  g. 

(b)  e  lies  entirely  in  H^UHf.  In  this  case  e  is  visible  from  g. 

(c)  e  crosses  one  or  both  planes  11^,17,  .  In  this  case  e  is  split  into  at  most 
three  subsegments,  one  of  which  {eg  =  er\H^nHf)  is  obscure  to  g,  and 
the  other  one  or  two  (e  —  eDH^nHf)  are  visible  from  g. 

So  far  we  have  assumed  that  the  faces  of  B  are  transparent,  but 
obviously  they  are  not.  However  the  effect  of  the  faces /^/y  incident  to  e  on 
the  visibility  between  e  and  g  is  completely  symmetrical  to  the  effect  of  the 
faces  /t  //  incident  to  g .  Therefore  they  can  induce  one  obscure  interval  on  g 
(gnHfnHj),  so  that  outside  this  interval  the  faces  of  B  do  not  affect  the 
visibility  of  e  from  g.  Hence  we  conclude  that  there  is  one  fixed  (possibly 
empty)  interval  eg  along  e  obscure  to  any  point  in  g-gHMfdHJ,  which 
proves  our  lemma.  □ 

The  above  lemmas  enable  us  to  calculate  in  time  O(n^)  all  pairs  (et,ej)  of 
edges  of  A  and  of  B  respectively  which  are  not  totally  obscure  to  one 
another,  and  for  each  such  pair  compute  the  two  corresponding  subsegments 
of  obscurity  obsi  C  e,  and  obsj  C  ej  as  in  Lemma  3.5.  (K  e,,ey  are  fully 
visible  from  one  another  then  obst  =  obsj=0.  )  Let  VS  denote  the  set  of  all 
such  quadruples  <ei,ej,obsi,obSj>. 

Nevertheless,  since  we  intend  to  use  VS  to  search  for  the  take-off  and 
landing  points,  we  must  also  take  into  account  the  possibility  that  the  take-off 
point  coincides  with  the  source  point  sq  itself,  or  that  the  landing  point 
coincides  with  the  destination  point  Z.  We  therefore  modify  VS,  by  inserting 
into  it  tuples  representing  visibility  of  edges  of  B  from  sq,  and  tuples 
representing  visibility  of  edges  of  A  from  Z.  Tuples  of  the  first  kind' have  the 
form  <so,ej,0,obsj>  where  obsj  is  the  obscurity  subinterval  of  ej  from  sq, 
provided  by  Lemma  3.4;  tuples  of  the  second  kind  have  the  corresponding 
symmetric  form  <ei,Z,obst,{)>. 

Note  that  if  the  points  t  and  Z  are  visible  from  one  another  then  /  =  Z, 
and  that  if  the  points  /  and  sq  are  visible  from  one  another  then  t  =  sq. 
Hence  we  can  delete  from  VS  all  tuples  <et,ej,obsi,obsj>  for  which  obsi 
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contains  the  obscure  subinterval  of  c/  with  respect  to  Z,  or  symmetrically  obsj 
contains  the  obscure  subinterval  of  ej  with  respect  to  sq.  Since  calculation  of 
the  obscure  subinterval  of  an  edge  with  respect  to  a  point  or  to  another  edge 
can  be  done  in  constant  time,  it  follows  that  VS  can  be  calculated  in  0{rr) 
time  (using  0{n^)  storage). 

Hence,  after  VS  has  been  computed,  the  implementation  of  the  crude 
algorithm  proposed  at  the  end  of  section  3.3  is  straightforward:  For  every 
tuple  <ei,ej,obsi,obsj>  (henceforth  e,  and  ej  may  denote  sq  and  Z  as  well, 
unless  otherwise  stated)  we  determine  the  length  of  the  shortest  path  which 
takes-off  from  a  point  in  the  interval (s)  ei  —  obsi  and  lands  on  a  point  in  the 
interval(s)  ej  —  obsj.  The  solution  of  the  general  problem  corresponds  to  that 
tuple  which  minimizes  this  path  length,  among  all  0{n^)  tuples  in  VS. 

Note  that  the  complexity  of  the  whole  procedure  depends  upon  the  time 
complexity  of  a  single  step  during  which  we  calculate  the  shortest  path  for  a 
fixed  tuple  (i.e.  for  a  fixed  pair  of  edges  of  A  and  B).  If  either  e,  or  ej 
degenerates  to  a  single  point  (^o  or  Z  respectively),  then  this  step  can  be 
performed  in  0{n)  time.  However  for  cases  where  ei,ej  are  both  edges,  each 
divided  into  0{n)  Voronoi  subintervals,  our  crude  algorithm  requires  0{n^) 
time,  to  perform  an  exhaustive  search  through  all  pairs  of  these  subintervals. 
Our  next  goal  is  to  improve  this  complexity  bound,  trying  to  reduce  it  as 
close  as  possible  to  0{n  log  n).  To  accomplish  this,  we  make  use  of  a  new 
kind  of  3-D  Voronoi  diagram,  which  we  call  Peeper's  Voronoi  diagram,  to  be 
introduced  and  analyzed  in  the  next  section. 

4.   Shortest  Path  Calculation  Using  Peeper's  Voronoi  Diagrams 

4.1.  Peeper's  Voronoi  Diagram 

Let  <e,e'  ,obs,obs'>  denote  a  fixed  tuple  in  VS  for  which  we  want  to 
determine  the  shortest  path  taking-off  from  a  point  in  e  —  obs  and  landing  on 
a  point  in  e'  —obs' .  For  simplicity  of  presentation  we  will  assume  in  this 
section  that  obs=obs'  =0,  so  that  the  edges  e,e'  are  entirely  visible  from 
each  other.  (K  obs, obs'  are  nonempty,  then  we  have  to  repeat  the  procedure 
described  below  at  most  four  times,  separately  for  each  pair  of  a  subsegment 
of  e  —  obs  and  a  subsegment  of  e'  —obs' .) 

The  edge  e  is  divided  into  0{m)  Voronoi  intervals,  induced  by  the  slices 
on  A.  Let  p^,  .  .  .  ,p^  denote  the  source  images  corresponding  to  these 
slices.  Similarly  the  edge  e'  is  divided  into  0{n)  intervals  induced  by  the 
slices  on  B,  with  q\,  .  .  .  ,q„  being  the  corresponding  images  of  the 
destination  point  Z.  We  claim  that  an  exhaustive  checking  of  all  0{n})  pairs 
of  intervals  on  e  and  e'  as  potential  locations  for  the  take-off  and  landing 
points  is  unnecessary,  because  most  of  these  pairs  will  not  induce  a  geodesic 
path  between  them.  This  fact  is  implied  by  the  following  argument.  Let 
Vor{qi,e')  denote  the  z-th  Voronoi  interval  along  e' .  Then  the  intervals  along 
e  which  can  match  Vor{qi,e')  in  the  above  sense  are  only  those  whose  source 
images  lie  nearest  (i.e.  yield  a  shortest  path  taking-off  from  e  )  to  some  point 
inside  Vor(qi,e')  . 
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Hence  we  wish  to  partition  e'  into  a  new  kind  of  Voronoi  subintervals, 
formed  by  a  certain  modified  3-D  Voronoi  diagram  of  the  set  of  points 
Pi,  .  .  .  ,p^  .  This  diagram  is  defined  so  that  if  a  point  x^e'  falls  inside  the 
cell  of  point  pi  in  the  diagram,  then  the  shortest  path  from  sq  constrained  to 
pass  through  e  \.o  x,  takes-off  from  the  i-th  Voronoi  interval  along  e.  The 
Voronoi  diagram  we  are  looking  for  differs  from  the  standard  3-D  diagram, 
because  the  paths  connecting  the  source  images  to  x  are  not  straight,  but 
rather  must  all  bend  over  the  edge  e.  Each  such  path  (being  geodesic)  must 
also  subtend  equal  incoming  and  outgoing  angles  at  its  point  of  contact  with 
e. 

Note  that  every  such  path  leaving  e  (at  a  point  t)  travels  in  space  along  a 
straight  line  until  it  reaches  e'  (at  a  point  /).  The  lines  e  and  //  define  a  half 
plane  P"^  (which  of  course  depends  on  /),  bounded  by  the  line  containing  e. 
Let  us  fix  such  a  half-plane  P  \  assume  without  loss  of  generality  that  e  lies 
on  the  Y-axis,  and  rotate  P"^  about  e  until  it  coincides  with  the  plane  IT 
containing  the  source  images  p^^,  .  .  .  ,p„  so  that  they  all  lie  in  the  left  half- 
plane,  while  P^  coincides  with  the  right  half-plane.  In  this  setting  all  the 
shortest  paths  taking-off  from  e  within  P"^  become  straight  lines  crossing  e 
from  the  left  to  the  right. 

Definition:  The  2-D  Peeper's  Voronoi  diagram  PV(S,e)  of  a  set  of  points  S 
(lying  in  the  left  half  plane)  and  an  edge  e  (lying  along  the  Y-axis),  is  a 
partition  of  the  right  half  plane  P"^  into  Voronoi  cells,  such  that,  for  each 
point  a^S,  the  corresponding  cell  V{a)  is  the  set  of  all  points  x^P'^  for 
which: 

(1)  ax  crosses  e  . 

(2)  |cu:|<|^a;|  for  every  b^S  for  which  bx  crosses  e  . 

In  addition,  P"*"  may  also  contain  a  cell,  denoted  as  V(0),  of  points  x  which 
do  not  belong  to  any  of  the  cells  V{a),a  ^S. 
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Figure  4.1  A  Peeper's  Voronoi  Diagram 

The  cells  of  PV{S,e)  may  have  rather  irregular  shape  and  in  fact  need  not 
even  be  connected  (see  e.g.  cells  V{b),V{c)  in  Fig.  4.1).  However  the 
boundary  between  different  cells  of  the  diagram  still  consists  of  finitely  many 
line  segments,  each  of  which  is  either  (part  of)  a  perpendicular  bisector  of  a 
segment  connecting  two  points  in  5,  or  else  is  (part  of)  a  ray  emerging  from 
an  end  point  of  e  and  lying  on  a  line  that  passes  through  a  point  in  S. 

Revolving  PV{S,e)  about  the  y-axis,  we  obtain  the  revolved  peeper's 
Voronoi  diagram  of  S  with  respect  to  e,  denoted  as  RPV{S,e).  It  consists  of 
conic  surfaces  of  revolution,  whose  generators  are  the  (straight)  edges  of  the 
2-D  peeper's  diagram. 

Returning  to  the  concrete  case  at  hand,  we  wish  to  calculate  the 
intersection  of  RPV{{pi,  .  .  .  ,/7^},e)  with  e' ,  thereby  obtaining  a  partition  of 
e'  into  new  peeper's  Voronoi  intervals  (not  to  be  confused  with  the  former 
Voronoi  intervals  Vor{qi,e'),  .  .  .  ,Vor{q„,e')  ).  For  each  such  subinterval  /, 
all  the  shortest  paths  from  sq  taking-off  from  e  to  the  points  in  /  are  induced 
by  the  same  source  image.  Our  goal  is  to  estimate  the  number  of  such 
subintervals,  and  to  provide  an  efficient  procedure  for  their  calculation.  For 
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this,  we  need  the  following: 

Lemma  4.1:  Let  S={a,b}  consist  of  just  two  points.  Then  the  edges  of 
PV{S,e)  lie  on  just  5  lines. 

Proof:  Let  the  endpoints  of  e  be  a:  1^2  •  The  five  lines  containing  all  edges  of 
PV{S,e)  are  the  four  lines  containing  the  segments  axi,  0x2,  bx^,  bxi, 
respectively,  emd  the  perpendicular  bisector  of  ab.  It  is  clear  that  the  position 
of  xiP"^  with  respect  to  each  of  these  lines  fully  determines  the  peeper's 
Voronoi  cell  {V{a),V{b)  or  V(0))  containing  x,  hence  each  Voronoi  edge  in 
PV{S,e)  must  lie  in  one  of  these  lines.  □ 

Remark:  A  closer  analysis  of  PV{S,e)  shows  that  it  consists  of  at  most  seven 
edges  (see  Fig.  4.2),  but  we  shall  not  make  use  of  this  property. 
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Figure  4.2  Peeper's  Voronoi  Diagram  of  two  Points 


Remark:  The  irregular  shape  of  the  peeper's  Voronoi  diagram  makes  it  more 
difficult  to  estimate  its  complexity.  It  is  an  open  problem  whether  the 
complexity  of  that  diagram  is  linear  in  the  number  n  of  points  in  5 ,  or  is  even 
close  to  linear.  K  this  were  the  case,  then  the  number  of  peeper's  Voronoi 
subintervals  on  e'  would  also  be  linear  (or  close  to  linear)  in  n,  and  the 
complicated  analysis  which  now  follows  would  then  become  redundant. 

Next  we  review  a  few  basic  facts  about  Davenport-Schinzel  sequences, 
which  we  will  need  to  obtain  a  sharp  bound  on  the  number  of  peeper's 
Voronoi  subintervals  induced  on  e'  by  e. 

Definition:  Let  n,s  be  positive  integers.  A  sequence  U={ui,  •  •  •  ,u„)  of 
symbols  is  an  (n,s)  Davenport-Schinzel  sequence  (a  DS(n,s)  sequence  for 
short),  if  it  satisfies  the  following  conditions: 

(1)   Each  Ui  is  one  of  n  allowed  different  symbols. 
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(2)  For  each  i<m  we  have  Ui¥'Ui+i  . 

(3)  There  do  not  exist  s  +  2  indices  l<ii<z2<  •  •  •  <is+2—'^  such  that 
Ui=Ui=Ui=  ■  ■  ■  =a  ,  ui=ui=ui=  •  •  •  =Z)  ,  and  a-^b  . 

Sharir  has  recently  obtained  an  almost  linear  upper  bound  on  the  length 
Xj(n)  of  {n,s)  Davenport-Schinzel_  sequences  (see  [Sh85a],[HS84]).  This 
bound  is  of  the  form  0(na(n)'^^'*("^'  ^)  ,  where  a(n)  is  the  functional  inverse 
of  Ackermann's  function,  and  is  thus  extremely  slowly  growing.  These  upper 
bounds  improve  a  previous  bound,  obtained  by  Szemeredi  [Sz74],  of  the  form 
X^(n)  <  C,n  log*/!,  where  log*Ai  is  the  height  of  the  smallest  exponential 

tower  2^      exceeding  n . 

Theorem  4.2:  The  number  of  intersections  of  e'  with  the  surfaces  of  the 
revolved  peeper's  Voronoi  diagram  RPV{{pi,  .  .  .  .Pm}>e)  of  e  and  its  source 
images  is  at  most  Xio(^  +  l)  • 

Proof:  Let  U=(ui,  •  •  •  ,ui^)  denote  the  sequence  in  which  u,  is  the  source 
image  whose  cell  in  RPV{{pi,  .  .  .  ,Pm)y^)  contains  the  i-th  interval  in 
e'  C\RPV{{pi,  .  .  .  ,p„},e).  The  sequence  U  satisfies  the  following  conditions: 

(1)  u,6{pi,  .  .  .  ,P;„}U{0}  for  each  i. 

(2)  For  each  i<k  we  have  Ui¥'Ui+i  (since  two  adjacent  intervals  on  e'  do 
belong  to  two  different  cells  of  the  revolved  diagram;  we  ignore  here 
degenerate  cases  in  which  one  of  the  surfaces  in  RPV({pi,  .  ,  .  .p„},e)  is 
tangent  to  e' ,  since  in  these  cases  the  point  of  tangency  can  simply  be 
discarded). 

(3)  There  do  not  exist  12  indices  l</i<Z2<...</i2^^  such  that 
Ut=Ui^=Ui^=...  =  a  ,  Ui^=Ui=Ui=...  =  b  ,  and  ai^b.  Indeed,  if  this  were 
the  case,  then  the  intersection  of  e'  with  RPV({a,b},e)  would  have  to 
consist  also  of  at  least  12  intervals.  This  however  is  impossible,  because 
RPV{{a,b},e)  consists  of  portions  of  just  five  conical  surfaces  (generated 
by  the  five  lines  containing  the  edges  of  PV{{a,b},e)),  each  of  which 
intersects  e'  in  at  most  two  points,  thus  yielding  at  most  11  intervals  in 
e'nRPV{{a,b},e)- 

Hence  U  is  aD5'(m  +  l,10)  sequence  so  that  its  length  is  at  most  \iQ(m  +  l).  □ 

4.2.   Constructing  the  Peeper's  Voronoi  Intervals  Induced  on  e'  by  e 

Theorem  4.2  will  next  be  used  to  obtain  an  efficient  procedure  for  the 
calculation  of  these  peeper's  Voronoi  diagram  intervals,  using  a  Divide  & 
Conquer  approach.  Let  S  denote  the  set  of  source  images  for  e.  We  divide  it 
into  two  sets  Sx,Si  of  equal  size,  and  calculate  recursively  the  intersection  of 
e'  with  /?PV(>Si,e)  and  of  e'  with  /^V'(52,e),  obtaining  two  partitions  111,172 
of  e'  into  peeper's  Voronoi  intervals.  Then  we  merge  ir^  and  it2  forming  a 
refined  partition  it  of  e' . 

Let  /  be  an  interval  of  it  which  is  contained  in  an  interval  /^  of  iri  and  in 
an  interval  l-^  of  172-  Let  p  ,q^S  be  the  source  images  corresponding  to  /^  and 
I2  respectively  (it  is  also  possible  that  p  or  5  is  0).  Our  merging  procedure 
has  to  replace  /  in  it  by  a  collection  of  0(1)  subintervals,  such  that  the 
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shortest  paths  from  sq  through  e  to  all  points  x  in  such  a  subinterval  are  all 
determined  by  the  same  source  image  (p  or  q).  To  do  this,  we  simply 
calculate  RPV{{p,q},e),  which  consists  of  at  most  five  conical  surfaces,  and 
intersect  these  surfaces  with  /.  Repeating  this  procedure  for  each  I^ir,  we 
obtain  a  refined  partition  -n*  of  e'  which  is  clearly  the  required  intersection 
of  e'  with  RPV{S,e).  (Note  that  if  //'  are  two  adjacent  intervals  in  it,  it  is 
possible  that  after  refining  them  into  subintervals  in  the  manner  just 
described,  an  extreme  subinterval  of  /  and  its  adjacent  extreme  subinterval  of 
/'  are  both  induced  by  the  same  source  image;  in  this  case  we  simply  merge 
these  two  subintervals  into  a  single  interval.) 

The  time  complexity  T(m)  of  this  algorithm  applied  to  a  set  S  of  m 
source  images,  can  be  estimated  by  the  recursive  formula 

T{m)  =  2T{m/2)+M{m) 

where  M{m)  is  the  time  required  to  merge  the  two  subpartitions  iTi,'ir2  of  e' . 
But  clearly  this  time  is  proportional  to  the  number  of  intervals  in  both 
partitions,  which,  by  Theorem  4.2  is  O (m a (m)*^ ("('")  ))  .  Solving  this 
recursive  formula  we  obtain 

r(m)  =  C>(ma(m)^("("')')logm) 

4.3.   The  Improved  Algorithm 

At  this  point  we  have  obtained  two  different  partitions  711,772  of  the  edge 
e' ,  where 

TTi  consists  of  the  peeper's  Voronoi  intervals  induced  on  e'  by  e,  and 

TT2  consists  of  the  (regular)  Voronoi  intervals  induced  on  e'  by  the  slices 
of  polyhedron  B  (with  Z  as  a  source). 

Let  us  now  merge  tti  with  1T2  to  obtain  a  refined  partition  it,  and  let  /  be  an 
interval  of  it  which  is  contained  in  an  interval  /^  of  tti  and  in  an  interval  1 2  of 
172.  Let  p,q  be  the  source  and  destination  images  corresponding  to  /^  and  1 2 
respectively.  Then,  by  the  methods  described  in  subsection  3.2,  one  can 
determine  in  constant  time  whether  there  exists  a  geodesic  path  taking-off 
from  the  Voronoi  interval  of  p  on  e  and  landing  on  I2,  and  also  calculate  the 
length  of  such  a  path  if  it  exists. 

171,172  contain  (9(ma(m)'^(°'('") ))  and  0{n)  intervals  respectively,  so  that 
17  contains  at  most  0(na(n)^^°^"^ ))  intervals.  Hence,  in  conclusion,  we  can 
determine  the  shortest  path  taking-off  from  e  and  landing  on  e'  in  the  time 
r(n),  which  is  almost  proportional  to  0{n  log  n),  and  thus  almost  achieve  the 
goal  we  set  at  the  end  of  Section  3.  Combining  this  with  the  overall  approach 
presented  in  section  3,  we  obtain  our  main  result: 

Theorem  4.3:  Given  two  convex  polyhedra  A,B  with  at  most  n  vertices  on 
each,  the  shortest  path  from  a  source  point  on  polyhedron  A  to  a  destination 
point  on  polyhedron  B,  which  does  not  penetrate  into  either  polyhedron,  can 
be  determined  in  0{n^a.{ny^'^^"^\ogn)  time  (and  0{n^)  storage). 
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5.   More  General  Problems 

In  sections  3  and  4  we  have  analyzed  shortest  paths  where  the  source  and 
destination  points  lie  on  the  surface  of  two  obstacles.  In  this  section  we  will 
consider  possible  extensions  of  this  case,  by  allowing  the  source  and 
destination  points  to  be  positioned  anywhere  in  free  space  surrounding  the 
two  convex  polyhedral  obstacles.  We  conclude  by  commenting  on  the  major 
difficulties  arising  when  trying  to  extend  our  technique  to  the  case  involving 
3  or  more  convex  obstacles. 

5.1.   Arbitrary  Source  and  Destination  Points 

Consider  two  disjoint  convex  polyhedra  A  ,B  in  3-space.  Let  sq  denote  the 
source  point  lying  on  the  surface  of  polyhedron  A  and  suppose  that  the 
destination  point  Z  is  positioned  in  free  space  surrounding  AJB.  Let 
Hull{A,Z)  denote  the  convex  hull  of  A  and  Z.  This  convex  hull  can  be 
computed  in  0{n)  time,  applying  a  merge  procedure  for  two  nonintersecting 
convex  hulls  (see  e.g.  [PH77]).  We  can  also  test  in  0{n  log  n)  time  whether 
the  intersection  of  the  two  convex  polyhedra  Hull{AX)  and  B  is  empty  (see 
[MP78]). 

If  Hull(A^)r\B  is  empty  then  clearly  the  shortest  path  tt  from  sq  to  Z 
(which  does  not  penetrate  into  A  or  B)  does  not  depend  on  B.  Hence  we  can 
ignore  B  and  consider  A  as  a  single  obstacle  in  space.  If  sq  and  Z  are  visible 
from  each  other  (which  can  now  be  determined  in  0(1)  time)  then  tt  is 
simply  the  straight  segment  sqZ.  Otherwise  we  replace  polyhedron  A  by 
Hull(A,Z)  and  use  Mount's  algorithm  for  finding  the  shortest  path  along  its 
surface  from  ^o  to  Z. 

If  Hull(A^)nB  is  not  empty,  we  use  the  following  Lemma: 

Lemma  5.1:  Let  Z  be  an  exterior  point  to  AyB  and  suppose  Hull{A^)nBi=0. 
Then  the  required  shortest  path  it  is  the  shortest  path  from  the  source  point 
Sq  (on  A)  to  the  point  Z  (on  the  convex  polyhedron  Hull{BX))  which  avoids 
penetration  into  either  of  these  two  polyhedra. 

Proof:  The  claim  of  the  lemma  is  based  on  the  observation  that  if 
Hull{A;Z)f\B^0  then  Hull{B^)r\A=0.  To  see  this  let  X^//u//(A,Z)n5. 
Then  there  exists  Y^A  such  that  X  lies  in  the  straight  segment  YZ.  Since  A 
and  B  are  disjoint  and  convex,  they  can  be  separated  by  some  plane  P.  But 
then  P  must  intersect  the  segment  XY  and  thus  must  also  separate  A  and  Z. 
Hence  P  separates  A  and  Hull{BX),  as  asserted. 

The  lemma  will  now  follow  immediately  if  we  show  that  the  shortest 
path  IT  never  enters  the  interior  of  Hull(B,Z).  There  are  two  cases  to  be 
considered: 

(1)  TT  passes  along  the  surface  of  B.  Let  r  denote  the  point  on  tt  where  the 
path  takes-off  from  B  towards  Z.  Then  r  must  lie  on  the  surface  of 
Hull{BX),  or  else  taking-off  towards  Z  could  have  been  accomplished 
earlier  along  tt,  thereby  making  it  shorter. 

(2)  TT  does  not  touch  B.  Let  t  denote  the  point  on  tt  where  the  path  takes-off 
from  A  directly  towards  Z.  Hence  points  t  and  Z  are  visible  from  each 
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other.  But  then  it  is  easy  to  see  that  tZ  does  not  meet  the  interior  of 
Hull{B^),  because  t  itself  lies  outside  this  hull,  and  because  tZ  does  not 
intersect  B .  □ 

Lemma  5.1  reduces  the  problem  of  finding  the  shortest  path  for  an 
arbitrary  point  in  space,  to  the  problem  studied  above  of  finding  a  shortest 
path  to  a  destination  point  lying  on  the  surface  of  another  polyhedron.  Next 
suppose  that  both  sq  and  Z  lie  in  the  free  space  surrounding  A,B.  We 
distinguish  between  two  cases: 

(a)  sq,Z  are  visible  from  each  other.  Then  'tt=sqZ. 

(b)  sq,Z  are  obscure  to  each  other.  For  this  case  we  have  the  following 
Lemma: 

Lemma  5.2:  Let  sq,Z  be  two  points  external  to  A,B,  which  are  not  visible 
from  each  other.  Then  the  shortest  path  it  from  sq  to  Z  which  avoids 
penetration  into  A  ,5  is  either  the  shortest  path  from  the  point  sq  lying  on  the 
surface  of  the  convex  polyhedron  HuII{A,sq)  to  point  Z  in  free  space,  or  the 
shortest  path  from  sq  lying  on  HuII{B,sq)  to  point  Z  in  free  space. 

Proof:  By  the  observation  made  in  the  proof  of  Lemma  5.1,  either 
Hull(A,so)nB  =  0  or  Hull{B ,SQ)nA  =  0.  Without  loss  of  generality,  suppose 
Hull{A,so)nB=0.  Then  Z  cannot  lie  in  Hull{A,so),  or  else  sq  and  Z  would  be 
visible  from  one  another.  It  now  follows,  as  in  the  proof  of  Lemma  5.1,  that 
IT  cannot  enter  the  interior  of  HuII{A,sq),  for  otherwise  it  could  be  shortcut 
from  Sq  to  the  first  point  at  which  it  leaves  that  hull.  This  proves  our  claim.  □ 

Lemma  5.2  reduces  the  case  of  an  arbitrary  source  point  in  free  space  to 
the  previous  problem  discussed  in  this  subsection.  Hence  we  have  in 
summary 

Theorem  5.3:  Given  two  convex  polyhedra  A,B  with  at  most  n  vertices  on 
each,  the  shortest  path  between  any  two  points  lying  either  in  free  space,  or 
on  the  surfaces  of  these  polyhedra.  which  does  not  penetrate  into  A  and  B, 
can  be  determined  in  C>(/i^a(«)^("("^^)log/i)  time  (and  0{n^)  storage). 

5.2.  Three  Convex  Obstacles 

Let  A,B,C  denote  three  disjoint  convex  polyhedra,  let  sq  denote  a  source 
point  lying  on  the  surface  of  A,  and  let  Z  denote  a  destination  point  lying  on 
the  surface  of  C.  Let  tt  denote  the  shortest  path  from  sq  to  Z  avoiding  the 
interiors  of  AJ3,C,  and  constrained  to  pass  through  the  surface  of  B  (we 
disregard  here  the  case  where  it  does  not  touch  B,  since  this  case  has  been 
studied  before,  and  can  be  checked  separately). 

Qearly  ir  is  the  concatenation  of  the  following  subpaths: 


(a)  A  shortest  subpath  along  the  surface  of  A  from  sq  to  a  first  take-off 
point  t^ . 

(b)  a  straight  Hne  segment  (in  empty  space)  connecting  t^  to  a  first  landing 
point  Ig  on  the  surface  of  B . 
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(c)  a  shortest  subpath  along  the  surface  of  B  from  /^  to  a  second  take-off 
point  tg . 

(d)  a  straight  line  segment  (in  empty  space)  connecting  tg   to  a  second 
landing  point  /^  on  the  surface  of  C. 

(e)  a  shortest  subpath  along  the  surface  of  C  from  Ic  to  Z. 

As  in  the  case  of  two  obstacles,  the  points  t^JsyhJc  ^^^^Y  determine  the 
required  shortest  path.  However,  finding  these  four  crucial  points  seems  this 
time  to  be  more  difficult. 

Suppose  we  know  the  two  edges  on  A  and  C ,  containing  the  points  r^  ,1^ 
respectively.  Assume  also  that  we  have  determined  the  slices  in  which  t^  and 
Ic  are  contained.  Then  we  can  fully  determine  the  structure  of  the  first  and 
the  last  subpaths  of  -ir.  However,  the  source  and  destination  points  of  the 
subpath  of  TT  along  B  are  not  known  in  advance,  so  that  we  have  to  consider 
separately  all  possible  pairs  of  edges  of  5,  in  which  the  points  lB,tg  might  lie. 
Furthermore,  even  for  a  fixed  such  pair  e,  e'  of  edges,  we  cannot  use 
Mount's  technique  to  partition  B  into  pre-calculated  slices  and  to  use  them  in 
the  calculation  of  the  middle  portion  of  tt. 

As  a  matter  of  fact,  the  problem  we  face  here  can  be  formulated  as  that 
of  calculating  all  possible  sequences  of  edges  of  B  which  can  be  crossed  by 
shortest  paths  from  any  point  on  e  to  any  point  on  e' .  K  we  had  all  these 
sequences  available,  then  we  could  unfold  each  such  sequence  4  in  turn,  so  as 
to  make  the  middle  portion  of  the  corresponding  candidate  path  into  a 
straight  segment,  and  then  find  the  shortest  path  from  sq  to  Z  constrained  to 
contain  such  a  segment,  a  task  which  can  be  accomplished  in  much  the  same 
way  as  was  done  in  Sections  3  and  4. 

This  problem  of  calculating  all  possible  sequences  of  edges  crossed  by 
shortest  paths  along  the  surface  of  a  convex  polyhedron  is  studied  in  a 
forthcoming  paper  by  Sharir  [Sh85b],  where  it  is  shown  that  the  number  of 
such  sequences  is  polynomial  in  the  number  of  faces  of  the  given  polyhedron, 
leading  to  polynomial-time  algorithms  for  finding  shortest  paths  between  two 
points  amidst  a  collection  of  a  fixed  number  of  convex  polyhedral  obstacles. 
These  algorithms  extend  the  observations  made  above,  and  are  also  described 
in  [Sh85b]. 
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